function [A, B ] = generateHyperbolic(Nx, Ny, coeff_x, coeff_y, coeff_c)
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here

%% Elliptic Equation

cc =ones(Nx+3,Ny+3);
ww =zeros(Nx+3,Ny+3);
ss = ww;
nn=  ww;
ee = ww;

%% The Inner Points
for ii=2:Nx+1
    for jj= 2:Ny+1
        cc(jj,ii) = coeff_c;
        ww(jj,ii) = -coeff_x;
        ss(jj,ii) = - coeff_y;
        ee(jj,ii) = coeff_x;
        nn(jj,ii) = coeff_y;        
    end
end
%% North and South Boundary Points

for ii = 2:Nx+1
    %   coeff = nn(ii,Ny);
    cc(Ny+1,ii) = cc(Ny+1,ii) + nn (Ny+1,ii);
    nn (Ny+1,ii) =0;
    cc (2,ii) = cc(2,ii) + ss(2,ii);
    ss(2,ii) =0;
end

%% East and West Boundary Points
 
for jj = 2:Ny+1
    %   coeff = nn(ii,Ny);
    cc(jj,Nx+1) = cc(jj,Nx+1) + ee (jj,Nx+1);
    ee (jj,Nx+1) =0;
    cc(jj,2) = cc(jj,2) + ww (jj,2);
    ww(jj,2) =0; 
end

A = generatespy_index(Nx,Ny, cc,ss,ww,ee,nn);


%[ ccc,css,cww,cee,cnn ] = rrb_index(Nx,Ny,cc , ss, ww);
[ ccc,css,cww,cee,cnn,m_cne,m_cnw,m_cse, m_csw  ] = rrb_hyper(Nx,Ny,cc , ss, ww);

 B = generatespy9_index(Nx,Ny, ccc,css,cww,cee,cnn,m_cne,m_cnw,m_cse, m_csw );

end

